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Introduction 


This  is  the  final  report  on  our  work  under  ARO-SDI  Contract  DAAL03-90- 
G-0108,  May  1,  1990  through  April  30,  1993. 

The  major  results  are  in  two  areas: 

1.  Issues  related  to  designing  configurable  and  fault-tolerant  processor  ar¬ 
rays  such  that  even  if  some  processors  are  faulty,  a  fault  free  array 
can  be  constructed  using  the  healthy  processors.  Such  studies  are  mo¬ 
tivated  by  applications  of  Wafer  Scale  Integration  (WSI)  where  for 
example,  a  large  number  of  processors,  configured  in  the  form  of  a 
square  grid,  can  be  put  on  a  single  wafer.  Due  to  finite  yield,  some 
of  the  processors  are  likely  to  be  faulty.  In  such  a  case,  one  can  work 
around  the  faulty  processors  and  reconfigure  the  remaining  processors 
to  form  a  healthy  grid.  Thus,  reconfiguration  methodologies  can  be 
viewed  as  tools  to  increase  the  effective  yield  of  the  WSI  processing 
arrays.  The  general  models  that  we  have  explored  consist  of  a  set  of 
identical  processors  embedded  in  a  flexible  interconnection  structure 
that  is  configured  in  the  form  of  a  rectangular  grid.  In  particular,  dur¬ 
ing  this  reporting  period  we  studied  models  that  use  multiple  tracks 
along  every  grid  line  and  developed  efficient  algorithms  for  reconfig¬ 
uration  in  such  models.  We  have  also  developed  efficient  scheduling 
algorithms  for  such  processor  arrays  in  the  presence  of  communication 
delays.  Our  results  on  this  topic  are  summarized  in  Section  2. 

2.  Aspects  of  antenna  array  processing  relating  to  methods  for  dealing 
with  unknown  or  partially  known  noise  field  and  the  array  mani¬ 
fold.  Our  approach  to  array  processing  problems  such  as  directions- 
of-arrival  (DOA)  estimation  has  been  based  on  the  very  successful 
Subspace  Method.  These  methods  however  need  explicit  knowledge  of 
the  noise  field  and  the  array  manifold.  However,  in  many  practical 
situations,  this  information  may  be  unknown  or  only  partially  known 
and  we  have  investigated  several  approaches  to  overcoming  such  lim¬ 
itations.  In  one  approach,  we  used  an  instrumental  variable  method 
where  the  sample  to  sample  temporal  correlation  of  the  signal  is  much 
larger  than  for  the  noise.  Specifically,  for  a  large  enough  differential 
delay,  the  covariance  of  the  array  outputs  can  show  a  finite  signal 
component  but  has  a  zero  contribution  from  noise.  Thus  the  unknown 
noise  covariance  problem  has  been  side  stepped.  We  have  provided  a 
detailed  analysis  of  this  approach.  Another  area  of  invetigation  has 
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been  the  analysis  of  combined  effects  of  finite  samples  and  model  errors 
on  DOA  estimation.  This  led  to  a  doubly  weighted  subspace  fitting 
approach  where  column  weighting  was  used  to  optimize  finite  sample 
effects  and  row  weighting  for  model  errors.  We  also  investigated  tech¬ 
niques  for  incorporating  prior  knowledge  of  array  response  errors  in 
the  estimator  using  a  Bayesian  approach  and  an  auto  calibraion  pro¬ 
cedure.  The  technique  here  is  to  allow  the  array  response  vectors  to 
adapt  to  given  data  within  the  limits  of  the  known  error  bounds  on  the 
array  manifold.  Our  results  on  this  topic  are  summarized  in  Section 
3. 

Two  Ph.D.  theses  in  these  areas  were  completed  during  this  contract 
and  several  papers  published,  accepted  and  submitted  -  see  pp.  13. 


2  Summary  of  Our  Work  on  Reconfigurable  Ar¬ 
rays 

Many  computations  in  matrix  algebra  can  be  conveniently  carried  out  on 
an  array  of  identical  processing  elements.  WSI  technology  provides  an  inex¬ 
pensive  approach  to  building  such  arrays.  However,  during  the  fabrication 
process  or  during  operation,  some  of  the  processing  elements  in  a  large  array 
are  inevitably  going  to  be  faulty.  Spare  PEs  and  extra  routing  hardware  are 
often  provided  so  that  a  fault-free  array  can  be  constructed;  such  reconfig¬ 
uration  capability  can  be  used  to  increase  the  yield,  and  to  guarantee  fault 
tolerance  in  applications  when  failure  is  not  permissible.  Our  work  in  this 
regard  has  been  concerned  with  the  design  and  analysis  of  such  configurable 
fault- tolerant  arrays  and  we  have  built  on  previous  work  within  our  group 
in  this  area,  fault-tolerant  arrays  and  we  have  built  on  previous  work  within 
our  group  in  this  area. 

The  general  model  considered  by  us  is  (see  e.g.  ,  [LSS89,  M0086,  Ros83, 
SS86,  SS82,  Sny82]):  it  consists  of  a  set  of  identical  processors  embedded 
in  a  flexible  interconnection  structure  that  is  configured  in  the  form  of  a 
rectangular  grid.  Each  grid  line  in  the  mesh  has  a  fixed  number  of  data 
paths  that  can  be  routed  along  it  (i.e.  ,  the  model  has  fixed  channel  width); 
switches  can  be  placed  at  every  grid  point  and  at  every  location  where  a 
processor  is  connected  to  the  grid.  Furthermore,  often  the  processors  are 
divided  into  a  set  of  non-spare  PEs  (say  an  m  x  n  array)  and  a  set  of  spare 
PEs  that  are  distributed  in  a  pre-determined  fashion. 

Given  a  set  of  faulty  PEs,  the  objective  is  to  reconfigure  the  connections 
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among  the  PEs  such  that  a  new  rectangular  logical  array  is  formed  compris¬ 
ing  only  the  healthy  PEs  and  demanding  no  more  hardware  resources  (e.g., 
spare  PEs,  tracks,  and  switches)  than  available.  It  is  clear  that  the  more  the 
additional  hardware,  the  higher  is  the  reconfiguration  probability.  Neverthe¬ 
less,  space  and  cost  limitations  place  practical  limitations  on  the  amount  of 
spare  PEs,  tracks  amd  switches  we  can  have.  A  question  that  now  arises  is: 
Given  a  configurable  architecture  with  fixed  resources,  are  there  efficient  and 
simple  algorithms  for  reconfiguring  such  architectures  with  high  probability? 

It  is  also  easy  to  see  that  if  the  number  of  faulty  PEs  is  less  than  the 
number  of  spare  PEs,  then  one  can  always  define  a  set  of  compensation 
paths  for  successful  reconfiguration.  A  compensation  path  is  defined  by  the 
sequence  of  substitutions  made  to  replace  a  faulty  PE  by  a  healthy  one. 
However,  the  characteristics  of  the  compensation  paths  (e.g.,  the  geomet¬ 
rical  distances  between  consecutive  nodes,  or  the  relative  positions  of  the 
nodes  in  the  grid)  determine  the  amount  of  routing  hardware  needed  to  im¬ 
plement  the  necessary  connections  among  the  logical  neighbors.  It  can  be 
easily  shown  that  if  the  number  of  routing  tracks  is  fixed,  then  one  cannot 
allow  arbitrary  sets  of  compensation  paths.  In  other  words,  by  limiting  the 
hardware  resources  one  limits  the  number  of  faulty  patterns  that  one  can 
reconfigure.  Hence,  another  natural  question  is  how  many  tracks  should  one 
provide  so  as  to  allow  a  large  enough  class  of  compensation  paths,  and  yet 
keep  the  hardware  redundancy  low. 

Yet  another  important  question  is:  given  a  network  of  healthy  (reconfig¬ 
ured)  processors  and  a  task  dependency  graph  with  communication  delays, 
what  should  be  the  schedule  of  the  tasks  so  as  to  maximize  performance. 

We  summarize  our  answers  to  these  questons/problems. 

2.1  New  Results 

We  developed  reconfiguration  techniques  for  a  3-track- 1-spare  model:  this 
consists  of  a  single  row/column  of  spare  PEs  on  the  side  of  the  array.  There 
are  three  routing  tracks  running  along  the  grid  lines  and  multiple  pro¬ 
grammable  switches  placed  at  the  intersections  of  the  tracks.  The  faulty 
PEs  are  not  assumed  to  be  able  to  function  as  interconnection  elements. 

We  developed  a  set  of  sufficient  conditions  for  reconfiguration  in  our 
model.  More  precisely  we  proved  that  if  there  exists  a  set  of  compensation 
paths  subjected  only  to  the  constraints  of  continuity  and  nonintersection, 
then  the  3-track- 1 -spare  model  can  always  accommodate  such  a  set.  We 
must  note  here  that  the  above  result  provides  the  first  known  theoretical 
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justification  of  the  observations  made  by  several  researchers  about  the  power 
of  3-track  models  (see  [JLS88,  YS89]). 

We  also  developed  polynomial  lime  algorithms  for  determining  such  a  set 
of  continuous  and  non  intersecting  compensation  paths  (if  it  is  possible  to  do 
so)  for  a  given  array  with  faulty  PEs.  If  there  is  no  such  set  of  compensation 
of  paths  that  covers  all  the  faulty  PEs,  then  the  algorithm  returns  a  set  of 
paths  to  cover  the  maximum  number  of  faulty  PEs.  Extensive  simulations 
were  carried  out  for  the  performance  of  our  3  track-l-spare  model  and  we 
compared  it  to  other  reconfiguration  models  that  use  similar  amount  and 
kind  of  additional  hardware.  The  results  show  that  our  model  performs 
better  than  other  models  that  use  the  same  or  greater  hardware. 

We  also  developed  efficient  algorithms  for  reconfiguration  in  a  neigh¬ 
borhood  reconfiguration  model  (see  also  [CF90,  SS86]).  In  this  model,  the 
neighborhood  reconfiguration  consists  of  an  n  x  k  rectangular  array  with 
one  column  of  spare  PEs  on  one  side.  The  routing  hardware  used  for  re¬ 
configuration  could  either  be  a  multiplexer  based  interconnection  scheme 
or  the  routing  is  implemented  by  tracks  and  switches.  The  goal  is  to  re¬ 
configure  the  array  minimizing  at  the  same  time  the  geometric  distances 
between  logical  neighbors  ( i.e PEs  that  are  connected  in  the  reconfigured 
array).  This  criterion  is  motivated  by  the  fact  that  shorter  interconnects 
reduce  the  communication  delays  among  the  PEs,  and  also  might  lead  to 
less  routing  hardware.  The  algorithm  that  we  present  has  100%  reconfigu¬ 
ration  yield  and  performs  provably  better  that  other  algorithms  that  use  the 
same  model  (see  for  example  [CF90,  FR85]).  We  developed  reconfiguration 
algorithms  for  the  case  of  only  one  column  of  spa-re  PEs  along  one  side  of 
the  array.  The  algorithms  are  based  on  a  simple  procedure  to  reconfigure 
an  jY  x  (A'  +  1)  array  into  an  (AT  -I-  1)  x  N  array  and  vice  versa.  We  also 
derived  time  complexity  of  these  new  algorithms. 

Finally  we  have  new  results  on  scheduling  problems  when  there  are  com¬ 
munication  delays  in  the  routing  network.  The  model  that  we  consider 
consists  of  m  identical  processors  P\ ,  P?, . . . ,  Pm  where  m  is  a  parameter  of 
the  problem  and  bounded  by  a  constant  c.  There  exist  n  computational 
tasks  Ti,T2,  . .  .,Tn  that  can  be  executed  in  any  of  the  m  available  proces¬ 
sors.  There  is  a  partial  order  among  the  tasks.  This  implies  if  Tj  — *  Tj  (i.e. 
Tj  is  dependent  on  Tj)  then  the  computation  of  Tj  in  some  processor  P:  can¬ 
not  start  before  task  Tj  is  computed  (say  in  processor  Pi)  and  the  results  of 
this  computation  have  been  transmitted  from  processor  Pj  to  processor  Pj. 
The  partial  order  among  the  tasks  introduces  a  precedence  graph  associated 
with  the  scheduling  problem.  The  tasks  {Tj}  are  the  nodes  of  the  graph  and 
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we  assume  directed  edges  between  the  nodes  X)  — *  T:  whenever  there  is  a 
dependency  between  tasks  T,  and  Tj  the  way  it  is  described  above.  We  shall 
assume  that  the  precedence  graph  is  of  the  out-forest  (or  in-forest)  form. 
All  the  tasks  are  of  unit  computational  time,  t.e.  it  takes  one  unit  of  time 
to  compute  any  of  the  tasks  T\ , . . . ,  Tn  in  any  of  the  processors  Pi,. .  .,Pm. 
The  processors  are  fully  connected,  t.e.  any  processor  can  communicate 
with  any  other  processor.  Moreover,  the  computation  of  the  tasks  in  the 
processors  is  independent  of  the  communication  among  them.  This  implies 
that  the  processors  can  be  running  tasks  at  the  same  time  that  communica¬ 
tion  is  taking  place  among  them.  The  communication  delay  among  any  two 
processors  takes  unit  time. 

For  such  problems  we  have  proved  the  following  theorem: 

Theorem  Given  out-forest  G,  it  can  be  determined  in  0(n2m)  time 
whether  there  exists  a  schedule  for  G  of  certain  given  length  A.  If  such 
a  schedule  exists  then  we  can  determine  this  schedule  in  0(n2m)  time. 

The  analysis  of  In-forest  precedence  graph  scheduling  problem  reduces 
to  the  Out-forest  precedence  graph  scheduling  problem  by  just  inverting 
the  direction  of  the  dependencies  between  the  tasks  and  then  inverting  the 
resulting  optimal  schedule.  In  the  process  of  proving  this  theorem  we  de¬ 
veloped  efficient  ways  of  transforming  the  given  graphs  that  are  subject  to 
communication  delays  into  delay  free  graphs  that  can  be  scheduled  with¬ 
out  taking  into  consideration  communication  delays  between  the  processors, 
and  whose  optimal  schedules  obey  the  precedence  and  communication  delay 
constraints  of  the  original  gTaph,  and  have  the  same  length  as  the  optimal 
schedule  of  the  original  graph.  This  allowed  us  to  use  dynamic  programming 
to  obtain  polynomial  time  algorithms  for  computing  the  optimal  schedule. 
Our  algorithm  is  a  generalization  of  the  results  presented  in  [DW85]  for 
graphs  with  no  communication  delays. 

Since  the  above  work  has  resulted  in  one  of  the  very  few  rigorous  results 
on  scheduling  with  communication  delays,  it  has  received  considerable  at¬ 
tention  in  the  research  community.  For  example,  Prof.  Eugene  Lawler  (of 
University  of  California  at  Berkeley)  [Law93]  has  used  some  of  the  concepts 
first  introduced  in  our  work  to  design  efficient  polynomial  time  algorithms 
for  determining  approximate  schedules;  our  algorithms  on  the  other  hand 
derive  optimal  schedules. 
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3  Summary  of  our  Work  ir  Array  Processing 

3.1  Introduction 

In  many  practical  signal  processing  problems,  the  objective  is  to  estimate 
from  noisy  measurements  a  set  of  constant  parameters  upon  which  the  under¬ 
lying  true  signals  depend.  For  example,  estimating  the  directions-of-arrival 
(DOAs)  of  impinging  wavefronts  given  the  set  of  signals  received  at  an  an¬ 
tenna  array  is  important  in  fields  such  as  radar,  sonar,  electronic  surveillance 
and  seismic  exploration.  High  resolution  frequency  estimation  is  important 
in  numerous  applications  including  Doppler  radar  and  system  identification. 
The  quantities  to  be  estimated  are  parameters  ( e.g .,  DOAs  of  plane  waves, 
nsoid  frequencies)  upon  which  the  observations  depend,  and  these  parame¬ 
ters  are  assumed  to  be  constant  over  the  observation  interval  which  is  long 
enough  to  collect  sufficient  data  to  ensure  parameter  estimates  of  the  desired 
accuracy. 

There  have  been  several  approaches  to  such  problems  including  the  so- 
called  Capon’s  maximum  likelihood  (ML)  method  [Cap69]  and  Burg’s  max¬ 
imum  entropy  (ME)  method  [Bur75].  Though  often  successful  and  widely 
used,  these  methods  have  certain  fundamental  limitations  (esp.  bias  and 
sensitivity  in  parameter  estimates),  largely  because  they  use  an  incorrect 
model  (e.g.,  AR  rather  than  special  ARMA)  of  the  measurements.  Pisarenko 
[Pis73]  was  one  of  the  first  to  exploit  the  structure  of  the  data  model,  doing 
so  in  the  context  of  estimation  of  parameters  of  cisoids  in  additive  noise 
using  a  covariance  approach.  Schmidt  [Sch79]  and  independently  Bienvenu 
and  Kopp  [BK79]  were  the  first  to  do  this  in  the  case  of  sensor  arrays  of 
arbitrary  form.  Schmidt,  in  particular,  accomplished  this  by  first  deriving  a 
complete  geometric  solution  in  the  absence  of  noise,  then  cleverly  extending 
the  geometric  concepts  to  obtain  a  reasonable  approximate  solution  in  the 
presence  of  noise.  The  resulting  algorithm  was  called  MUSIC  (Multiple 
Signal  Classification)  and  has  been  widely  studied  since  its  inception.  The 
geometric  concepts  upon  which  MUSIC  is  founded  form  the  basis  for  a  much 
broader  class  of  subspace-based  algorithms  and  is  the  focus  of  this  discussion. 

Since  the  pioneering  work  of  Schmidt,  several  new  techniques  based  on 
the  subspace  approach  have  been  developed.  Notably,  the  ESPRIT  tech¬ 
nique  (Estimation  of  Signal  Parameters  via  Rotational  Invariance  Tech¬ 
niques)  proposed  by  Paulraj,  Roy  and  Kailath  [PRK85],  [PRK86j,  [RK89]. 
More  recently,  new  results  in  multi-dimensional  techniques  such  as  the  clas¬ 
sical  maximum  likelihood  method  and  Weighted-Subspace-Fitting  (WSF)  of 
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Ottersten  and  Kailath  [VOK89]  have  attracted  attention  due  to  their  po¬ 
tentially  superior  perfomance.  In  fact,  an  optimal  choice  of  weighting  can 
be  shown  to  yield  estimates  that  achieve  the  Cramer- Rao  lower  bound  on 
the  error  variance  [OVK92].  Techniques  such  as  MUSIC,  ESPRIT  and  WSF 
are  collectively  known  as  subspace  methods.  These  methods  assume  that  the 
noise  field  is  spatially  white  and  completely  known  to  within  a  scale  factor. 
Likewise,  the  response  of  the  array  to  a  wavefront  from  any  given  angle  (also 
known  as  the  Array  Manifold)  is  also  assumed  to  be  known. 

3.2  New  Problem  Areas 

However  there  are  in  practice  several  situations  where  the  above  assumptions 
may  not  be  true.  Due  to  changes  in  weather,  the  surrounding  environment, 
and  antenna  location,  the  response  of  the  array  may  be  significantly  different 
than  when  it  was  last  calibrated.  Furthermore,  the  calibration  measurements 
themselves  are  subject  to  gain  and  phase  errors.  For  the  case  of  analytically 
calibrated  arrays  of  nominally  identical,  identically  oriented  elements,  er¬ 
rors  result  since  the  elements  are  not  really  identical  and  their  locations  are 
not  precisely  known.  Depending  on  the  degree  to  which  the  actual  antenna 
response  differs  from  its  nominal  value,  the  algorithm  performance  may  be 
significantly  degraded.  The  requirement  of  known  noise  statistics  is  also 
difficult  to  satisfy  in  practice,  since  the  surrounding  environment  and  orien¬ 
tation  of  the  array  may  be  time-varying.  In  addition,  one  is  often  unable  to 
account  for  the  effect  of  unmodeled  “noise”  phenomena  such  as  distributed 
sources,  reverberation,  noise  due  to  the  antenna  platform,  and  undesirable 
channel  crosstalk.  Measurement  of  the  noise  statistics  is  complicated  by  the 
fact  that  there  are  often  signals-of-interest  observed  along  with  the  noise 
and  interference.  Because  of  these  difficulties,  it  is  often  assumed  by  de¬ 
fault  that  the  noise  field  is  isotropic,  that  it  is  independent  from  channel 
to  channel,  and  that  its  power  in  each  channel  is  equal.  When  the  SNR  is 
high,  deviations  of  the  noise  from  these  assumptions  are  not  critical  since 
they  contribute  little  to  the  array  covariance.  However,  at  low  SNR,  the 
degradation  may  be  severe. 

We  briefly  discuss  below  our  results  in  handling  cases  where  the  above- 
mentioned  assumptions  of  known  array  response  and  noise  statistics  is  not 
correct.  Earlier  work  has  been  previously  reported  in  [Swi91],  [Fri90],  [SH91] 
Several  authors  have  also  investigated  DOA  estimation  algorithms  that  at¬ 
tempt  to  mitigate  the  effects  of  the  antenna  and  noise  model  errors  described 
above  [PK85,  PK86,  PRSK86,  FW88,  WWL88,  WF88,  W0V91], 
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Main  Results 

A  variety  of  models  could  be  used  to  describe  the  perturbed  array  matrix. 
In  practice,  the  response  of  a  given  sensor  is  typically  known  to  within  some 
tolerance  in  gain  and  phase  that  accounts  for  variations  in  the  construction 
of  the  sensor  and  the  conditions  under  which  it  is  to  operate.  This  tolerance 
may  be  specified  as  limits  above  and  below  some  nominal  response,  or  as  an 
expected  deviation  around  the  nominal.  Consequently,  one  might  assume 
that  A  is  specified  in  probabilistic  terms  ( e.g the  mean  and  variance  of 
the  elements  of  A  are  assumed  known),  and  that  the  sensor  array  is  just 
one  realization  from  the  probability  space  of  arrays.  A  particularly  simple 
model  of  this  type  that  has  been  widely  used  [Fti90,  LV90,  ZW88,  Kur89] 
is  to  assume  that  the  array  rsponse  vectors  are  zero-mean,  white,  circu¬ 
lar,  and  stationary:  T.  the  errors  are  independent  from  sensor-to-sensor,  its 
covariance  is  clearly  diagonal.  Off-diagonal  terms  indicate  sensor-to-sensor 
correlations  that  result,  for  example,  if  there  are  uncalibrated  mutual  cou¬ 
pling  effects,  or  if  some  sensors  tend  to  perturb  uniformly  (such  as  identical 
or  adjacent  elements).  We  have  derived  (VS]  exact  expressions  for  the  DOA 
estimation  error  covariance  of  several  popular  subspace  based  methods  as¬ 
suming  the  error  model  given  above  in  the  presence  of  finite  data  samples. 
In  particular,  this  is  done  for  a  weighted  version  of  MUSIC  and  subspace 
fitting  algorithms  using  row  and  column  weightings.  An  important  outcome 
of  our  analysis  is  that,  for  a  given  error  model,  proper  weightings  lead  to 
minimum  variance  DOA  estimates. 

Accurate  signal  parameter  estimation  from  sensor  array  data  is  a  prob¬ 
lem  which  has  received  much  attention  in  the  last  decade.  A  number  of 
parametric  estimation  techniques  have  been  proposed  in  the  literature.  In 
general,  these  methods  require  knowledge  of  the  sensor-to-sensor  correla¬ 
tion  of  the  noise,  which  constitutes  a  significant  drawback.  This  difficulty 
can  be  overcome  only  by  introducing  alternative  assumptions  that  enable 
separating  the  signals  from  the  noise.  In  some  applications,  the  raw  sen¬ 
sor  outputs  can  be  pre-processed  so  that  the  emitter  signals  are  temporally 
correlated  with  correlation  length  longer  than  that  of  the  noise.  An  Instru¬ 
mental  Variable  (IV)  approach  can  then  be  used  for  estimating  the  signal 
parameters  v  hout  knowledge  of  the  spatial  color  of  the  of  the  noise.  We 
have  developed  [VS093]  a  improved  technique  that  can  give  significantly 
better  performancerformance. 

As  pointed  earlier,  the  subspace  methods  require  an  exact  character¬ 
ization  of  the  array,  including  knowledge  of  the  sensor  positions,  sensor 
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gain/phase  responses,  mutual  couplings,  and  receiver  equipment  effects.  Un¬ 
less  all  sensors  are  identical,  this  information  must  typically  be  obtained  by 
experimental  measurements  (calibration).  In  practice,  of  course,  all  such 
information  is  inevitably  subject  to  errors.  Recently,  several  different  meth¬ 
ods  have  been  proposed  for  alleviating  the  inherent  sensitivity  of  parametric 
methods  to  such  modeling  errors.  The  technique  proposed  herein  is  related 
to  the  class  of  so-called  auto-calibration  procedures,  but  it  is  assumed  that 
certain  prior  knowledge  of  the  array  response  errors  is  available.  This  is  a 
reasonable  assumption  in  most  applications,  and  it  allows  for  more  general 
perturbation  models  than  does  pure  auto-calibration.  The  optimal  maxi¬ 
mum  a  posteriori  (MAP)  estimator  for  this  problem  has  been  formulated, 
and  a  computationally  attractive  large-sample  approximation  derived.  We 
have  shown  [VS93]  that  our  approach  is  statistically  efficient,  and  verified 
this  through  computer  simulation. 
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